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p |. ABSTRACT 

O ' Aims. We present a new measurement of the cosmic X-ray background (CXRB) in the 1 .5-7 keV energy band, performed by exploiting 

,Q , the Swift X-ray telescope (XRT) data archive. We also present a CXRB spectral model in a wider energy band ( 1 .5-200 keV), obtained 

t/^ . by combining these data with the recently published Swift-BAT measurement. 

, ^ , ■ Methods. From the XRT archive we collect a complete sample of 126 high Galactic latitude gamma-ray burst (GRB) follow-up 

' observations. This provides a total exposure of 7.5 Ms and a sky-coverage of ~7 square degrees which represents a serendipitous 

, survey, well suited for a direct measurement of the CXRB in the 1.5-10 keV interval. Our work is based on a complete characterization 

^ ■ of the instrumental background and an accurate measurement of the stray-light contamination and vignetting calibration. 

' Results. We find that the CXRB spectrum in the 1.5-7 keV energy band can be equally well fitted by a single power-law with photon 

' index r=1.47±0.07 or a single power-law with photon index r= 1.41 ±0.06 and an exponential roll-off at 41 keV. The measured flux 

I in the 2-10 keV energy band is 2.18±0.13 xlO"" erg cm^^s^'deg-^ in the 2-10 keV band. Combining Swift-XRT with Swift-BAT 

I . (15-200 keV) we find that, in the 1.5-200 keV band, the CXRB spectrum can be well described by two smoothly-joined power laws 

• ' with the energy break at 29.0±0.5 keV corresponding to a vF,, peak located at 22.4±0.4 keV. 

'"^ ' Conclusions. Taking advantage of both the Swift high energy instruments (XRT and BAT), we produce an analytical description of 

, the CXRB spectrum over a wide (1.5-200 keV) energy band. This model is marginally consistent with the HEAOl measurement 

OO . (-10% higher) at energies higher than 20 keV, while it is significantly (30%) higher at low energies (2-10 keV). 
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$_i 1. Introduction 90% ) CXRB, the resolved fraction decrea sing at higher ener- 

9^. r^u ■ ^ ^ , J .^vr,T,^ • „ j j gies ("Mo rettl et al.ll2003l: lmrslev et al."2005). The i?oint sources 

The Cosmic X-ray background (CXRB) is usually defined Cueing the resolved fraction of CXRB in the 2-10 keV band 

as the integrated emission of all the extragalactic sources m ^^^^ ^^^^ ^^^^^ ^^^^j ^.^^ a small contributions 

the X-ray energy band (-2-100 keV). The name background from gal axy cluste rs and star burst galaxies (iHorn schemeier et d! 

comes directly froin the first X-ray astronomical observatioii 200O; Ba TieT^ 2004; Br andt & Hasin^eJ FoOS: TozzT^ 

(| GiaccometalJ | l962|) , when an apparendy diffuse background 20O6). Furthermore, a highly anisotropic diffuse component is 

^^1^^^' , extra-solar X-ray source ^^^^^^ .^^ j^^^^ ^j^^^^ ^ j^^^ (Sohan 200?.). This is con- 

(Sco X- 1 ). The CXRB spectral properties, flux and isotropy were ^.-^^^^^ ^ Local Bubble, the Gala ctic h"al5lGaleazzi et all 



accurately (10%) measured over a wide energy band by the A2 ^mjH and the intergalactic medium dciK&OittiEB [HH), 



and A4 experiments on board the High Energy Astronomical ^j^jj^ ^^ ^^^ .^^^ ^.^ q^I^^^.^ ^^^^^^ ^^^^ 

Observatory 1 (HEAOn_sgteUite. The analytical model pro- component is negligible, 
duced by IGruber et alj ( 119991) . combining A2 and A4 observa- 
tions with higher energy data has been considered as a refer- There is a general consensus on the sources from which the 

ence for many years (G99 model hereafter). Much effort has CXRB origina tes, and the bac kground paradox can be consid- 

been spent to quantify the fraction of CX RB emission due to un- ered solved dSetti & Woltjei<1989); nevertheless, the spectrum of 

resolv ed point sources. As predicted by i Giacconi & Zamoranil the X-ray integrated emission is still very important in the study 

(Il987h a combination of large and deep surveys performed of the statistical properties of those sources that are too faint 

by focusing telescopes in the soft part of the X-ray spectrum to be detected individually by currendy operating telescopes, 

(< 10) keV has succeeded in resolving almost the entire (80- as highly absorbed AGNs and very high red shift quasars. The 

extrapolation of the AGN observed spectra (unabsorbed and 

Send offprint requests to: alberto.moretti@brera.inaf.it Compton-thin) to the region of the CXRB peak (~30 keV) can 
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Table 1. A compilation of CXRB flux measurements in the soft 
energy band sorted by year of publication, compared with the 
G99 model. 



Instrument 


Flux 2-10 keV 


Reference 




[ 1 0" ' ' erg cm"- s" ' deg"^ ] 




HEAOl 


1.65±0.17 


|Gruberetal.Ul999») 


Rockets 


2.20±0.20 


McCammon et al. 1 1983) 


ASCA-SIS 


1.92±0.09 


Gendreau et al. 1 1995i) 


SAX-MECS 


2.35±0.10 


Vecchi et al. ( 1991) 


ASCA-GIS 


1.94±0.20 


Kushino et al. 1 200?) 


RXTE-PCA 


1.64±0.05* 


Revnivtsev et al. (2003) 


XMM-Newton 


2.24±0.16 


De Luca & Molendi (2004) 


HEAOl -A2 


1.66+0.08* 


Revnivtsev et al. (2005) 


Chandra 


2.19±0.26 


IHickox & Mai-kevitclj a006) 



(*)The original values are 1.91±0.06 and 1.96±0.10 respectively. We 
correct them, accord ing to the values reported in Table 3 of 
dRevnivtsev et al.|[2005 ^ to account for the cross-calibration with 
XMM-Newton. 
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account for ~ 75% of the peak (iGiUi et al.l2007l) . Compton-thick 
AGNs are thought to be responsible for the remaining fraction 
dTreister & Urrvil2005l; lOiUi et al.ll2007l; iBallantvne & PapovichI 
I2007h . Given that even the most recent AGN surveys in the 
hard band, >10 keV, can a dd only few p ercent to this number 
(ISazonov et allllOOTl 120081: iTueller et al.l l2008). the CXRB pro- 
vides a key boundary condition in the determination of the cen- 
sus of the heavily obscured AGNs. Moreover, an accurate mea- 
surement of the CXRB spectrum is an important observational 
constraint in the study of the very high redshift (z>6) AGNs 
which will remain u nresolved even with t he next generation 
of X- ray telescopes dSalvaterra et alj|2007t iRhook & Haehneltl 
120081) . FinaUy, a proper characterization of the CXRB spec- 
trum is also crucial to ensure proper background subtraction 
in the study of low surface brightness diffuse X-ray emission 
coming from the outskirts of clusters and g roups of galaxies 
(iGastaldeUo et alJl2007l: ISnowden et alJl2008l) . 

Measurements performed in the soft part of the CXRB spec- 
trum with different instruments (see Table [T]) yield a scatter 
which is much larger t han the one expect ed from standard can- 
dle flux measurements dKirsch et al.ll2005l) . meaning that the ob- 
served discrepancy cannot be entirely explained by the differ- 
ences in absolute calibrations of the individual instruments. The 
large scatter and the poor control on systemati c uncertainties in 
the CXRB measurem ents led some authors (iUeda et al.l 120031: 
iTreister & Um^l2005h to use the G99 model, re-normahzed by a 
factor of ~30%. The underlying (but not verified) assumption is 
that the G99 spectrum is correct in shape but with the normaliza- 
tion aff ected by some c alibra t ion p roblems of the HEAOl instru- 
ments. IWorslev et alJ d2005L 12006) use an even more artificial 
solution, combining the X MM-Newton CXRB measurement 
dPe Luca & Molendill2004l) up to 8 keV and the re-normalized 
G99 at higher energies. On the other hand, recently published 
measurements, performed by means of wide-field not-focused 
hard X-ray instruments (SAX-PDS, INTEGRAL-IBIS, Swift- 
BAT) yield results co nsistent (10% level) with the G99 model in 
the 20-50 keVraii ge dChurazov et al.ll2007HF^ntera et al.ll2007t 
lAjello et alll2008l), rev ersing the recent trend that prefers higher 
intensities "(IUedall2007h . 

Here we present a new measurement of the CXRB spec- 
trum in the 1.5-7 keV energy band, obtained by the analysis 
of the archival data of the X-ray te lescope (XRT) on board the 
Swift satellite dGehrels et alJl2004 . a mission dedicated to the 



Fig. 1. Upper panel: Sky distribution of the 126 observations in 
Equatorial coordinate, with the Galactic plane exclusion (dotted 
lines). The size of the points is proportional to the exposure and 
it is not representative of the observed field size. Middle panel: 
the cumulative survey sky coverage as a function of (logarithm 
of) exposure time. Lower panel: distribution of the (logarithm 
of) nominal exposure time. 

study of gamma-ray bursts (GRBs) and their afterglows. XRT 
uses a Wo her I mirror set, or iginally designed for the JET-X 
telescope dCitterio et all 11994 . to focus X-rays (0.2-10 keV) 
onto a XMM-Newton/EPIC MOS CCD detector dBurrows et alJ 
2005). GRBs are detected and locahz ed by the Burst Alert 
Telescope (BAT)( Barthelmy et aljllool) . in 15-300 keV energy 
band and followed-up at X-ray energies (0.3-10 keV) by the 
X-Ray Telescope. Following-up gamma-ray burst afterglows, 
the Swift-XRT obtains deep exposures on random positions of 
the sky, totally uncorrelated with already known bright X-ray 
sources, providing us with a simple and direct measurement of 
the CXRB spectrum. 

2. Work strategy 

For each energy E we can consider the signal Stot registered in 
a typical high Galactic latitude GRB afterglow follow-up obser- 
vation, as the sum of 4 factors. These are the GRB signal, the 
CXRB itself, which is the one we aim to measure, the elec- 
tronic and particle induced background (NXB or instrumental 
background) and the stray-light (SL), i.e. the contamination from 
sources outside the telescope field of view. 

Stot(E) = CXRB(E) X fvign(E) + NXB(E) + SL(E) + GRB(E).(1) 

The GRB afterglow signal can be easily eliminated by filter- 
ing data both in space and in time. The CXRB itself is con- 
tributed by bright resolved plus faint unresolved sources and 
it is affected by vignetting (fvign)- We measure the NXB using 
two different and independent datasets: a two- day observation 
performed with the focal plane camera assembly (FPCA) sun- 
shutter closed (SC) and the data collected in a region of the de- 
tector which is not exposed to the sky (NES). We evaluate the 
third element, the SL contamination using a series of off-axis ob- 
servations of bright sources. Given the high level of the CXRB 
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Fig. 2. XRT detector. The inner dotted circle indicates the detec- 
tor region we consider for the CXRB measurement. The larger 
dotted circle indicates the nominal field of view. The continuous 
circle shows the conservative definition of field of view we use to 
define the NES regions. The continuous box is used to define the 
NES border and show the regions contaminated by calibration 
source (CS) out of time events. 



isotropy (IShafer & Fabianlll983l;lR"evnivtsev et al.ll2008l) we can 
consider this factor as a fraction of the CXRB (SL = fsLxCXRB) 
in such a way that: 



CXRB(E) = 



Stot(E) - NXB(E) 

fvi<:„(E) -)- fsL(E) ■ 



(2) 



To measure the CXRB spectrum we perform stacked spectral 
analysis of a large sample made only of GRB follow-up obser- 
vations. The following four sections are devoted to fully describ- 
ing the technical details of our work, which allowed us to esti- 
mate the terms of Eq. |2] and their uncertainties. In particular, 
in Section [3] we describe the dataset and the reduction proce- 
dures; in Section|4]and |5]we present how we measure the NXB 
and SL contamination, respectively. In Section |6] we describe 
how we calculate the vignetting correction and its uncertainty. 
In Section |2] we present the spectral analysis procedure and its 
results, discussing them in Section [8] 

Throughout this paper, all errors are quoted at 90% confi- 
dence level, unless otherwise specified. The photon index is de- 
noted as F. 



3. Sample selection and data reduction 

During the first 45 months of operation (Jan 2005- July 2008), 
the Swift-XRT observed some 300 GRB afterglows with typi- 
cal exposure times of 70-100 ks during the ~ 10 days follow- 
ing the prompt event. We consider all the long (T90 > 2s) GRB 
follow-up observations with a nominal standard exposure longer 
than 10 ks and Galactic latitude \b\ > 20°. Because we find long 
term variations of the NXB level we consider only data after 
January 2006 in such a way that the NXB scatter remains lower 
than 10% (see Section |4|i. Similar variations (<3% per year) in 
the NXB level were observed in SAX-LECS-MECS and ASCA- 



GIS (iPerri & Giommil2000l:lKushino et al.l2002l) and were inter- 
preted as due to a gradual drop in the satellite altitude and/or to 
the cycle of solar activity. 

The sample consists of 126 GRB observations from January 
2006 to July 2008. For each observation we exclude from our 
analysis the data collected in the first day (the segment 0) in or- 
der to exclude the brightest phases of the afterglows. For each 
observation we consider only the central 200 pixel radius (7.9 
arcminutes) circle, excluding a 30 pixel radius (1.18 arcmin- 
utes), corresponding t o ~95% of the enc ircled energy fraction 
of a point-like source dMoretti et al.ll2005l) around the GRB po- 
sition. This corresponds to a nominal field of view (FOV) of 
0.054 square degrees. The real observed sky solid angle varies 
from observation to observation depending on the precise point- 
ing distribution of the observation. 

We reduce data using the standard software (HEADAS soft- 
ware, v6.4, CALDB version Dec07) and following the proce- 
dures in the instrument user guide Q. We replace the standard 
good time interval (GTI) definition, which is tuned for the obser- 
vations of bright point-like sources, by more restrictive filters. 
Due to the failure of the thermo-electric cooler power supply, 
the XRT CCD temperature is subject both to orbital (4°C in 5.9 
ks) and long term (15°C on a d ay time scale) varia tions, rang- 
ing between -70°C and -47°C jKennea et al.ll2005h . Dark cur- 
rent and hot pixels are highly temperature dependent and create 
high instrumental background in the low energy band (0.3-0.7 
keV) d uring observations performed at temp eratures higher than 
-52°C ("Pagan i et al.ll2007l:lMoretti et alj200 7). Moreover, due to 
the low orbit of the Swift satelUte, a typical target can be ob- 
served no more than 1-2 ks on a single orbit. Therefore the data 
from single object are split in different segments. Occasionally 
some reflected light from the Earth limb significantly increases 
the very low energy (< 0.5 keV) background at the end or at the 
beginning of an observation segment. To reduce these effects, 
we select intervals with CCD temperature <-55°C and elevation 
angle (i.e. the altitude of the observation direction on the Earth 
horizon) > 40°, instead of the standard -47°C and 30°, respec- 
tively. Moreover, we consider only data from observation seg- 
ments longer than 300 seconds and eliminate the first and the 
last 100 seconds of each orbital segment. After the complete 
time-filtering procedure, these procedures typically reduce the 
effective exposure time to 50% of the standard ones. The total 
nominal exposure time considered is 7.5 Ms with median value 
of 40 ks for single observations. The final exposure time distri- 
bution of the 126 observations, together with the sky-coverage 
is shown in Fig.[T] the surveyed area is 7 and 1.3 deg^ at 10 and 
100 ks respectively. 

4. Instrumental and particle induced signal 

To evaluate the instrumental and particle induced background 
(NXB) we use two different datasets. 

First, we use the two day observation performed between 
2007-09-04 18:50:00 UT and 2007-09-06 18:42:00 UT when 
the instrument Sun shutter (0.38 mm thickness stainless steel of 
grade 302) was closed due to an improper slew which brought 
the XRT to point ~ 15 degrees from the Sun. The instrument au- 
tomatically closed the shutter in front of the camera. For the next 
two days the usual XRT observations were performed, but with 
the shutter closed. We apply to these data the same reduction 
and filtering procedures that we apply to the sky data. The final 
exposure time for the instrumental background with the shutter 



http://heasarc.nasa.gov/docs/swift/analysis/ #documentation 
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Fig. 3. NXB count rate (per pixel) in the 1.5-7 keV energy band 
as measured by SC and NES datasets. The linear gradient in the 
CCD vertical direction is evident. 




Satellite Time [s since 2001-01-01] 



Fig. 4. The instrumental Ni Kq. line flux registered in the 56 ob- 
servations with exposure times longer than 40ks. The continuous 
line is the average, while the dotted line is the best linear fit. The 
total scatter of these measurements with respect of the average 
is at the level of 10%(1 cr), as shown in the right panel of the 
figure. 



closed is 67 ks with an average count rate of 1.92+0.05x10 
(7.71+0.3x10-**) counts s-'pixer' in the 0.3-10 (1.5-7.0) keV 
energy band. For the remainder of the paper we will refer to this 
dataset as shutter closed (SC) dataset. 

The second dataset is provided by the data collected in 
the regions of the detector which are not exposed to the sky 
(NES). These are four different regions (2507 pixels each) 
close to the CCD boundary and delimited by the FOV and 
the corner sources (Fig. |2]|. The FOV and corner source re- 
gion definitions are reported in the standard calibration file 
(CALDB) swxregion20QlQl(81v003.fits. In particular the 
nominal field of view of the telescope is the 300 pixel radius 
circle centered on the detector pixel (300,300). We conserva- 
tively adopted a wider definition of the FOV which is the 323 
pixel radius circle centered in the detector pixel (307,300). Then, 
we define the NES as the parts of the box centered in detec- 
tor pixel (307,300), width 436 and height 596 lying outside the 
conservative FOVQ- The signal registered in these regions has 
been telemetered since June 2008 when the telemetered detector 
area was increased to 600x600 pixels. We use all the available 



^ In simpler words, according to ds9 syntax, this corresponds to 
box(307, 308,436, 596,8)-circle(307, 300, 323) 




OFF- -AXIS angle [aromin] 

Fig. 5. Upper Panel: off'-axis observations of the Crab nebula, 
used to calibrate the SL contamination. Lower Panel: ratio be- 
tween SL and on-axis flux from the Crab at different ofF-axis 
angles. 



data present in the Swift-XRT archive of the photon counting 
(PC) observations between June- July 2008. This results in a to- 
tal exposure of 2.4 Ms. The uncertainty in the determination of 
NXB is one of the main sources of uncertainty of this measure- 
ment. For both these datasets we assume that the signal is con- 
tributed only by the particle induced and the pure instrumental 
background. Possible sources of systematic errors are the fol- 
lowing. First, the fluorescence from the shutter itself: if the shut- 
ter produces some fluorescence lines, our NXB estimate would 
be systematically higher than the correct value. Second, signif- 
icant inhomogeneities in the CCD response or in the intrinsic 
fluorescence background, could bias the NXB measurement in 
the NES regions. The third error source is the time dependence 
of the NXB. As already mentioned, in some previous missions, 
where the NXB background has been estimated by means of ob- 
servation of the dark Earth, a time dependence on time scale of 
years has been observed. 

To check our data for these systematic errors, we first verify 
consistency between the two datasets, SC and NES. Comparing 
the two datasets we find that the NXB in the 1.5-7 keV band 
displays a gradient in the vertical direction of the CCD, with 
the bottom regions being 30% lower than the top regions. This 
trend is very well described by a linear fit with SC and NES 
datasets being highly consistent (Fig.|3]l. The consistency of the 
two datasets give us good confidence that the first two sources 
of systematics are negligible for our purposes. We do not have a 
direct way to monitor the NXB time dependence during all the 
observations (NES data started to be telemetered only in June 
2008) over the entire energy band. Nevertheless we can use the 
data in the 7-8 keV interval, where the CXRB signal is low and 
the detected signal is dominated by the Ni K„line at 7.48 keV 
produced by the fluorescence of the telescope material. This is 
uniformly distributed over the detector area with a typical count 
rate of 3.8x10"** count s"' pixel"'. We compare the Ni line ob- 
served in the NES regions with the one observed during the SC 
and sky observations. To minimize the statistical error we con- 
sider the 56 observations longer than 40ks. We model the data 
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Fig. 6. Variations of the effective area as function of the off-axis 
angle, usually called vignetting. The dashed line represents the 
standard calibration function, while the black points are the val- 
ues we find for the 2-3 keV energy band. 



in the 6.8-8.2 keV energy band by means of a Gaussian plus a 
straight line. We find that the Ni line flux decreases slightly with 
mission time (correlation at 2.5cr confidence) producing scatter 
of ~ 10% with respect to the average (Fig. |4|. We note that the 
line flux registered in SC and NES data is perfectly consistent 
with the average of the observed fluxes. For the purpose of the 
stacked spectral analysis we account for this uncertainty in the 
final error budget as explained below in Sect. |2l 



5. Stray-light 

The other main component of the non-cosmic background is the 
SL. Ray-tracing simulations indicate that this is produced by 
photons coming from sources that are outside the telescope FOV 
at distance between 20 and 100 arcminutes from the optical axis 
of the telescope, whereas the FOV of the telescope (mirror + 
detector) has a radius of ~ 12 arcminutes. A fraction of the pho- 
tons produced by these sources reach the detector after only one 
reflection on the mirror or even directly, passing through the mir- 
ror shells without any interaction. Some X-ray telescopes mount 
baffles on top of the mirrors to prevent such a contamination . 
This is not the case for XRT, for which the SL is a significant 
fraction of the diffuse radiation registered on the CCD. 

In order to evaluate the level of contamination in XRT im- 
ages, we take advantage of the many Crab Nebula calibration 
observations. Then a series of observations at 7 different off-axis 
angles, ranging from 15 to 90 arcminutes were performed. We 
compare them with the on-axis calibration observation. We cal- 
culate the fraction of the source flux present on the central 200 
pixel radius circle of the detector as the ratio between the flux ob- 
served at each distance from the optical axis and the on-axis flux. 
Given the large off-axis angles the dimensions of the Crab neb- 
ula are negligible for our purposes. The results of this analysis 
are shown in Fig.|5] We split our analysis in several different en- 
ergy bands, finding no significant variation as function of energy 
up to 5 keV. These observations clearly show that XRT images 
are contaminated by the emission of sources outside the field of 
view up to ~ 70 arcminutes. Due to the isotropy of the CXRB we 
can calculate the expected contamination as the surface integral 
of the (measured) relative flux produced by the Crab observed at 
different off-axis angles. The result of the integration in the 1 .5-7 
keV band is the CXRB fraction fsL= 0.268+0.015 (see Eq.|2li. 




200 



400 600 
Energy Channel 
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Fig. 7. The different components of the signal registered during 
an observation with no bright source present in the field of view. 
The black line is the whole signal (Stot) registered during the 
35 ks observation of the afterglow of GRB080319C. Data have 
been reduced and filtered as explained in Section [3] The red line 
is the NXB component (NES data). The green line is the ex- 
pected SL contamination. In the lower panel the ratio between 
the total NXB (instrumental +particle induced + SL) and the Stot 
is plotted. The XRT energy channel are 0. 1 keV wide. 



To give an idea of the SL contam ination in absolute te rms, 
using the LogN-LogS calculated by iMorettiet al] ( l2003h . we 
find that for each XRT image we expect a contaminating flux 
of 2.7+0. Ix 10^'^erg s"' diffuse over the 200 pixel central cir- 
cle in the 1.5-7 keV band, corresponding to a count rate of 
6.3+0.2x10"^ (assuming a spectral photon index of 1.4). 



6. Vignetting 

Because the CXRB is an extended source with a uniform surface 
brightness profile, the variation of the effective area as function 
of the off-axis angle, i.e. the vignetting, must be accounted for. 
For each energy the vignetting correct ion can be analytically 
described by a polynomial function (iTagliaferri et al.l |2004|) . 
Therefore, first, we calculate the vignetting correction as func- 
tion of the off-axis angle using the standard calibration (CALDB 
coefliicients swxvign2001OlOlvO01 . fits). Then, the total vi- 
gnetting factor (/vign(E) in Eq. [TJ is given by the integration 
of this function over the 200 pixels radius circular region. In 
the 1.5-7 keV energy band the integrated vignetting correction 
ranges from 6% to 14%. 

To estimate the vignetting calibration accuracy, which is a 
relevant factor in the uncertainty calculation of our measure- 
ment, we consider all the point-like sources (1945 sources) de- 
tected in our survey, excluding GRB afterglows. We calculate 
the median count rate of these sources in different off-axis angle 
bins. Assuming that the source populations detected at differ- 
ent off-axis angles coincide (once we eliminate GRB afterglows) 
and having the necessary statistics to make cosmic variance neg- 
ligible, the ratio between the median count-rate on- and off-axis 
give us the vignetting coefficient. We repeat this operation for 
different energy bands and we find that the standard vignetting 
calibration is accurate at the level of a few percent (Fig.|6]l. 
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Fig. 8. CXRB spectrum and the power law best fit. Best fit values 
are reported in Table |2] 



7. Spectral analysis 

The XRT nominal energy band is from 0.3 to 10 keV. The frac- 
tion of Stot due to non-cosmic background (NXB-hSL) during a 
typical XRT observation (without any bright source in the FOV) 
depends on energy and is shown in Fig. [T] In the energy range 
between 1 and 6 keV the NXB contributes ~50% of the total 
signal registered during the observation. As a comparison, in the 
XMM-Newton data the cos mic fraction of the total signal is 20% 
dPe Luca & Molendill2004l) . Above 7 keV the NXB and, in par- 
ticular, the particle induced background becomes dominant, due 
to the presence of the Ni (K,, and Kp at 7.478 and 8.265 KeV) 
and Au (Lq, at 9.713 keV) fluorescence lines. We only consider 
the energy band between 1.5 and 7 keV, excluding the data with 
energy less than 1 .5 keV because the Galaxy contribution is not 
neghgible jKushino et al.ll2002l: Iffickox & Markevitchl2006l) . 

We account for the vignetting by mod- 
ifying the nominal ancillary response file 
(ARF,swxpc0tol2s0(6)2001ISlOlv010.ar£) by the /vign(E) 
factor (see Section |6]l . We modify the ARF file also to account 
for the SL contamination according to Eq. |2] with the results 
reported in Section |5] Finally, we calculate the contamination 
of the GRB afterglow residua ls outside the 30 p ixel radius 
using the analytical PSF model iMoretti et al.l (120051) . The GRB 
residual signal is energy dependent and contributes a maximum 
of 1.1% below 2 keV, becoming negligible above 3 keV. We 
account for this contamination applying another small energy 
dependent coiTection to the ARF file. To summarize, for each 
energy E we modify the nominal ARF according: 

ARF'(E) = ARF(E) x (fv.g„(E) + fsL(E)) x (1 - GRB,.es(E)). (3) 

In order to account for the CCD defects and the excluded 30 
pixel radius circular region around the GRB position we create 
an exposure map for each observation. The overall correction, 
weighted for the exposure time of the single observation, cor- 
respond to 6.3% for the sky exposure and 3.2% for the back- 
ground. We coiTect this by modifying the BACKSCAL keyword 
in the spectrum (PHA) files. For the background file we consider 
the NFS regions that provide better statistics than the SC ob- 
servations (Section HJ. The four NFS regions are not homoge- 
neous due to the spatial gradient in the NXB. Nevertheless the 
symmetry of the geometry allows us to use this dataset without 
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Fig. 9. Upper panel: the distribution of the CXRB flux measure- 
ments from the 113 observations with effective durations longer 
than 10 ks (black line), together with the best Gaussian fit (red 
line). The scatter in the measurements is the sum of the statisti- 
cal error (~15%) and cosmic variance. Lower panel: The com- 
parison between CXRB variance expected (continuous line and 
small circles) with the observed one, corrected for the statistical 
contribution (triangle). 



any coiTection. We perform the stacked spectral analysis merg- 
ing the 126 event files, re-sorting the events and the GTIs, and 
extracting the spectrum from the 200 pixel radius central circle 
in detector coordinates. To fit the data and calculate the fluxes 
we use XSPEC(vl2.4). 

To account for the systematic uncertainties in the NXB 
(Section |4]i, vignetting factor (fyign, Sect. |6|l and SL contam- 
ination (fsL, Sect. |5]l measurement we produce a large num- 
ber (10,000) of simulated datasets, randomly varying the NXB 
normalization, the fyign and fsL, according to the appropriate 
Gaussian distributions. For the NXB normalization we use the 
mean standard deviation that we observe for the Ni line fluxes in 
our sample (lower panel of Fig.|5]l. For the fyign and fsL we con- 
servatively use a standard deviation equal to 5% which slightly 
exceeds the estimated errors. 

We neglect the Galactic contribution (absorption and emis- 
sion) and fit our data by means of a simple power law, obtaining 
the numbers reported in Table |2] and shown in Fig. [H For anal- 
ogy with previous w orks in the literature tGruber et al J 1 1 999t 
iFrontera et al.ll2007l) . we also fit our data with a cut-off" power 
law (CPL) with the energy break fixed to 41.13 keV. Both the 
models provide a good description of our data in the 1.5-7 keV 
energy interval. 

Finally we combine our data with the Swift-BAT CXRB 
measurement perf ormed in the energy band 20-150 keV 
(Ai ello et al.l 12008) and we fit the joined energy distributions 
with two smoothly joined power laws (2S JPL) of the form 



dN 
dE 



C ■ E^ 



[keVcm^^ s"' deg"^] 



(4) 



This is the same model lAjello et alj ( l2008h uses to fit a large col- 
lection of CXRB measurements, together with the Swift-BAT 
new measurement. We note that, because a specific response 
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Fig. 10. CXRB measures plotted against time, Galactic latitude 
and CCD temperature. Here the plotted error are only the sta- 
tistical ones. We fit all the three datasets with a straight line 
(y=qx-i-a), finding that they are consistent with a constant (q=0) 
at the level of 1 cr. 



20081). Assuming the LogN-LogS calculated by iMoretti et"an 
(2003*), which is flatter at low energies and generating 1000 ran- 
dom samples with different dimensions (ranging from 0.01 to 1 
square degree of sky), in the flux range 10""'-10""' erg cm^^s"', 
we find that a more realistic value is Q"°^, as shown in the lower 
panel of Fig.|9] We find that, with 0.054 deg^ XRT field, we ex- 
pect a variance of 15.1%. This is very close to the one we ob- 
serve in our sample, once we account also for the contribution 
of the statistical error: cr^ = cr^ - cr^=14.3+1.8%. We note also 
that the extrapolation of our simulations to the total surveyed 
area (~7 deg^) tells us that the stacked analysis uncertainty due 
to cosmic variance is negligible. 

7.2. Data check 

We check our data against any bias due to the time of the obser- 
vation, CCD temperature, and Galactic latitude. As explained in 
Section |4] we find a slight dependence of the NXB on the time 
of the observations. However, as the NXB dependence is slight 
and the NXB contribution is minor the 113 CXRB flux measure- 
ments do not have any significant correlation with the observa- 
tion time (see lower panel of Fig. [TO] i. As already said (Section 
O the XRT CCD temperature is variable due to the fact that it 
is only passively controlled. Because dark current and hot pixels 
are temperature dependent, we also plot the 113 flux measure- 
ments against the average temperature of the observations (see 
central panel of Fig. [TO] i. Finally we check the flux measure- 
ments against Galactic latitude to exclude any significant contri- 
bution from our galaxy to the XRT measurements (upper panel 
of Fig. [Tol l. For all the three datasets we find that the best linear 
fit is consistent with a constant at the level of 1 cr . 



matrix has been produced in the Swift-BAT measurement, the 
cross-calibration factor for point-like sources jGodet et al.l2008l) 
cannot be applied here. 

Our best fit values are reported in Table |2] The peak of the 
energy distribution Epeak is given by 

Epeak = Eb-U—|1 =22.4 + 0.4 keV. (5) 
7.1. Cosmic variance 

To study the variance of our sample we consider the 113 obser- 
vations with durations longer than 10 ks for which the spectral 
parameters and the flux of the CXRB can be calculated with an 
acceptable accuracy (cts ~15%). We find that the flux distribu- 
tion is well described by a Gaussian with a standard deviation of 
(To=20.8%+2.4 (upper panel of Fig. |9]l. The maximum CXRB 
flux value in our sample, at ~ 5 ctq from the mean, is observed in 
the field of GRB 071028B where the galaxy cluster Abell SI 136 
is present. This is not surprising, given the fact that Abell clus- 
ters are ~5000 distributed over ~27,000 deg^ of sky meaning 
that 1-2 Abell clusters are expected in our survey. 

The variance we observe in the flux distribution is con- 
tributed by both by the statistical error (cr^) and cosmic variance 
((Tc). We find that the latter is consistent with the CXRB variance 
expected for the area surveyed by a single observation if we as- 
sume that the CXRB is entirely produced by point sources. In 
fact, it can be shown that, if we assume that the source fluxes are 
distributed as the classical F"^''^ LogN-LogS, the cosmic vari- 
ance scales with the surveyed area as Dr^-^, dRevnivtsev et al.l 



8. Discussion 

As mentioned in the Introduction, the CXRB spectrum normal- 
ization is still a debated issue. 

The Swift-XRT measurement, we present here, is very close 
to XMM-Newton (TablefTJ. This is not unexpected as the Swift- 
XRT effective area calibration has been slightly modified to 
match XMM-Newton by means of simultaneous observations 
0. For what concerns cross-calibration, Swift-XRT measures 
fluxes 5-10% lower t han RXTE-PCA d uring simultaneous ob- 
servations of 3c273 (iGodet et al.l l2008h . Cross-calibration ob- 
servations of IE 0102.2-7219, the brightest supernova remnant 
in the Small Magellanic Cloud, recently showed that Chandra- 
ACIS, XMM-Newton-MOS, Suzaku-XI S and Swift-XRT agre e 
to within +10% for all instruments dPlucinskv et all l2008l) . 
Therefore, the diff^erences with HEAOl and RXTE-PCA mea- 
surements cannot be entirely explained b y the absoulte calibra - 
tion differences, as already pointed out bv lFrontera et al.l (|2007|) . 

In the region of the CXRB peak (~30 keV) all the measure- 
ments show a acceptable agree ment(~10%), in the soft band 
the XRT measurement confirms iRevnivtsev et al.l (l2005l) con- 
clusions: in the 2-10 keV band narrow-field focusing telescopes 
measure CXRB values which are significantly higher than the 
ones found by wide-field not focusing telescopes. However, the 
XRT data, although inconsistently higher than the G99 model, 
smoothly join the higher energy data as we show by the good fit 
to the XRT and BAT. 

Below 60 ke V, the G99 model consists in a CPL with r= 1 .29 
and energy break 4 1 . 1 3 (note that no uncertainties are reported in 

3 SWIFT-XRT-CALDB-09-Vn available at 

|http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/index.html| 
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Table 2. The best fit results for the three diff'erent models we use. Among brackets are reported the statistical contribution to the 
total error budget. 
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Fig. 11. Upper panel: A compilation of flux measurements both in the soft and hard energy bands. F or the clearness of the plot, not 
all the soft energy measurements reported in Table [T] are shown here. Because 'Grub er et al.l (1 19991) do not report the uncertainties 
in the best fit values, we use a fiducial 5% error for G99 model. Lower Panel: Ratio of the flux measurements plotted in the upper 
panel with our joined XRT-hBA T fit. Colors are the same of upper panel. For comparison, with the dotted line we plot (only in the 
bottom panel) the ad-hoc model iWorsley et aP (l2005i) used to calculate the resolved fraction. 



the lGruberet al.l (1 19991) paper). As shown in Fig. [TT] the differ- 
ences from the G99 model range from 30% below 10 keV down 
to 5-10% in the region of the CXRB energy peak. This is due to 
the fact that the slope of the soft part in our best fit is steeper 
(1.41 instead of 1.29) and the peak of the spectrum is much 
softer (22 keV in stead of 29 keV). As previously discussed, 
lAiello et alj (2008) uses a 2SJPL to fit Swift-BAT data together 
with a large collection of different CXRB measurements down 
to 2 keV. In comparison to this model we find that the soft en- 



ergy slope is significantly softer (1.41+0.02 versus 1.32+0.02), 
while the high energy slope, the energy break and the normal- 
ization are consistent. Interestingly, our model has the same 
CPL shape, with energy break at 41.13 keV an d photon index 
1.4, t hat provides the best fit to SAX-PDS data (iFrontera et al.l 
I2007h in the 20-50 keV band, albeit with a significant difference 
in normalization. Finally we also note that XRT data and our 
model are well consi stent with the INTEGRAL measurement 
dChurazov et al ] l2007l) all over the considered energy band. 
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In summary, breaking the paradigm that G99 spectrum has 
the correct shape shows that CXRB data can be analytically de- 
scribed by a 2S JPL with the values reported Table|2]and a peak in 
the energy distribution at 22.9+0.4 keV. In the 1.5-50 keV range, 
this function is very similar to a CPL with the energy break of 
41.13 keV and photon index of 1.4. 

We note that the 2-10 keV CXRB flux measurement di- 
rectly affects the ca lculation of the CXRB resolved fraction. 
iMoretti et al.l (l2003b . combining shallow and deep surveys and 
integrating the source number counts, esti mate that the resolved 
fraction of 2-10 keV CXRB is 87+6% . IWorslev et al] ([2005h 
refined this calculation finding that the resolved fraction ranges 
from 80% in the 2-4 keV band to < 60% for energies higher 
than 6 keV. The main reason for the inconsistenc y between the 
two re sults is the value of the CXRB total flux. ' Moretti et alj 
( I2003h used an average of a sample of CXRB measurements, 
yielding a value of 2.02+0.11x10"" erg cm^^s^'deg"^ which 
is 10% less than the presen t measurement. As already men- 
tioned, IWorslev et al.l (l2005h used an ad hoc model, combin- 
ing the XMM-Newton measurement with a re-normalized G99 
model. As shown in the FigHT] (bottom panel) this model, al- 
though not motivated from an observational point of view, is not 
very far from our best fit. If we assume the present measurement 
for the CXRB and integrate the LogN-LogS of iMoretti et all 
(l2003l). we find a result for the CXRB resolved fraction which is 
79±6% in the 2-10 keV band, in very good ag reement with the 
average value quoted bv IWorslev et al.r(l2005l) . The values rela- 
tive to the single narrow bands at higher energies, on the other 
hand, should be slightly corrected, applying our CXRB value. 

The LogN-LogS extrapolation at very low fluxes (10"'^ 
erg cm^^s^'deg"^, a factor 20 lower than the faintest Chandra 
deep field sources) cannot account for all the CXRB. This im- 
phes that a not negligible fraction of th e CXRB is supposed 
to b e produced by non det e cted sources. IWorslev et all (|2006|) 
and iHickox & MarkevitchI (12007 ) correlate almost the entire 
CXRB unresolved fraction to optical/IR detected galaxies in 
the Chandra deep fields. These are star-forming galaxies which 
are expected to overwhelm the number of AGNs at very low 
fluxes (Ranaflietal. 2003; Bauer etal. 2004), absorbed AGN 
jTreister & Urrvll2005l: IGiUi et al.ll2007l) which are supposed to 
be the main component at higher energies and with a small con- 
tribution from very high redshift (z>6) quasars. 

9. Conclusion 

We use the Swift-XRT archival dataset to determine the flux 
and spectrum of the CXRB. This has two main advantages. 
The first one is the observational strategy which provides us 
with a truly random sampling of the X-ray sky, not correlated 
with previously known sources. The second is the low level of 
the NXB background, which allows measurement of the CXRB 
with high accuracy. Similar to other focusing telescopes, we find 
fliat CXRB flux is significandy higher than HEA01/G99 model. 
Nevertheless combining our dataset with Swift-BAT data, we 
show that we can describe the CXRB spectrum with a simple 
model (two smoothly joined power laws) over a wide energy 
band. The model we propose is much more observationally mo- 
tivated than the ones recently used in the literature for population 
synthesis models and for the CXRB resolved fraction calcula- 
tion. Using the present CXRB measurement we calculate that 
the resolved fraction in the 2-10 keV energy band is 79±6%. 
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